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Introduction 

Previous papers have used two important functions for the solution of fractional order 
differential equations, the Mittag-Leffler function E q [at q j (1903a, 1903b, 1905), and the 

F-function F q [a y t] of Hartley & Lorenzo (1998). These functions provided direct solution and 

important understanding for the fundamental linear fractional order differential equation and for 
the related initial value problem (Hartley and Lorenzo, 1999). 

This paper examines related functions and their Laplace transforms. Presented for 
consideration are two generalized functions, the R -function and the G -function, useful in 
analysis and as a basis for computation in the fractional calculus. The R -function is unique in 
that it contains all of the derivatives and integrals of the F-function. The R -function also returns 
itself on glh order differ-integration. An example application of the R -function is provided. A 
further generalization of the R -function, called the G -function brings in the effects of repeated 
and partially repeated fractional poles. 

Functions for the Fractional Calculus 

This section summarizes a number of functions that have been found useful in the solution 
of problems of the fractional calculus and more particularly in the solution of fractional 
differential equations. 


Mittag-Leffler Function 

The Mittag-Leffler (1903, 1903, 1905) function is given by the following equation 

£,[ ' ]= lf(^)' " >a 


0) 


This function will often appear with the argument -at q , its Laplace transform then, is given as 

L&i-ai'b i i « >0 - (2) 

1 ,l * larfa+OJ st?'+a) 
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AgarwaFs Function 

The Mittag-Leffler function is generalized by Agarwal (1953) as follows 



I 


m=0 



r {am + py 


( 3 ) 


This function is particularly interesting to the fractional order system theory' due to its Laplace 
transform, given by Agarwal as 

(9 

This function is the (a — /?) order fractional derivative of the F-function, (of Robotnov (1969) 
and Hartley (1998)), with argument a = 1, to be presented later. 


Erdelvi’s Function 

Erdelyi (1954) has studied the following related generalization of the Mittag-Leffler function 



~ m 

y £ 

{am + py 


a, (5 > 0, (5) 


where the powers of t are integer. The Laplace transform of this function is given by 


l{e.M= 


Y r(m + 1) 
Sr {am + py * 1 


a ,p >0. 


( 6 ) 


As this function cannot be easily generalized it will not be considered further. 


Robotnov and Hartley’s Function 

To effect the direct solution of the fundamental linear fractional order differential equation 
the following function was introduced (Hartley and Lorenzo, 1998) 





n~Q 


(~a)V« 
Y{nq + q ) 


q > 0. 


(7) 


This function had been studied earlier by Robotnov (1969, 1980) with respect to hereditary 
integrals for application to solid mechanics. The important feature of this function is the power 
and simplicity of its Laplace transform, namely 

’ 4 >0 - ( 8 ) 

s H -a 


Miller and Ross' Function 

Miller and Ross (1993, pp.80 and 309-351) introduce another function as the basis of the 
solution of the fractional order initial value problem. It is defined as the v th integral of the 
exponential function, that is 


E,(y,a) = 



= t' e a 'y' (v,at) = 


r 


Y— M_ 

sr(v+*+iy 


(9) 
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where y *(v, at) is the incomplete gamma function. The Laplace transform of equation (9) 
follows directly as 


L{£,(v,a)}= — — Re(v)> 1. 

S~Cl 


(10) 


Miller and Ross then show that 


4 x« H £,(/v-i,‘> 4 =-r— • ? = 1,2,3.., v = i = 

[P J s -a q 2 3 


CO 


which is a special case of the F-function of Robotnov and Hartley. 

The above functions are studied in considerable detail by their originators and others. The 
interested reader is directed to the supplied references. 

A Generalized Function 

It is of significant usefulness to develop a generalized function which when fractionally 
differintegrated (by any order) returns itself. Such a function would greatly ease the analysis 
of fractional order differential equations. To this end the following is proposed, consider 
the function 


R g . v [a,c,t] = ^ 


(a? (t-cf 


P r((n + l)?-v) ' 


( 12 ) 


Our interest in this function will normally be for the solution of fractional differential equations 
for the range of t > c = 0. For t <c, R will be complex except for the cases when the exponent 

((« + l )<7 - 1 - v) is integer. The more compact notation 


p r((n + l)?-vj 


(13) 


is also useful, particularly when c = 0. 

The Laplace transform of the R -function is 




WzfgLL = f (ax ' 1 

n=o r((n + i)?-v) „ =0 [r((n +i)?- v)J 


Consider first the case for c = 0, then we have 






Now from (Erdelyi et al, 1954) 

iff }=r(v + I>“ v '‘ Re(v)>— 1, Re(s)> 0. 

Then equation 1 5 becomes 

1 




n=0 




Re((n + l)^-v)>0, Re(s)>0. 


(14) 

(15) 

(16) 

(17) 
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Re((« + l)^-v)>0, Re(.s)>0. (l8) 



This can be written as a geometric series that converges when J a/ s q | < 1. It can be shown, by 
long division, that 

Re(<?-v)>0, Re(s)>0. (l9) 

Now for c & Othe shifting theorem for the Laplace transform (Wylie p. 281) is 

L{f(l-b),(l-b)}=e- l "L{f(l)} b> 0, (20) 


ife,(o.0,r)}= 


s q - a 


where the unit step function u{t - b) effectively causes /(/ - b) = 0 for t < b . Under the 
assumption that R q v [a,c,f]= 0 for t < c, this theorem and the result (equation 19) are applied 
to yield 

c > 0,Re((« + \)q - v)> 0,Re s > 0. (2l) 

Table 1, in a later section, presents a summary of the defining series and respective Laplace 
transforms for the functions discussed in this paper. 




^ ~ CS „ V 

e s 




Properties of the R q%v (a 9 c,t) Function 

The genera] time domain character of the R -function is shown in figures 1, 2, and 3. Figure 1 
shows the effect of variations in q with v = Oand a = ±1. The exponential character of the 
function is readily observed (see, q = 1 ). Figure 2 shows the effect of v on the behavior of the 
R -function. The effect of the characteristic time a is shown in figure 3. The characteristic time is 
M a q . For q = 1, 1 / a is the time constant, when q- 2 we have the natural frequency, when q 
takes on other values we have the generalized characteristic time (or generalized time constant). 




Figure la. Effect of q on R 0 (l,0,f), 
v = 0.0, <3—1.0 


Figure lb. Effect of q on R 0 (-l,0,r), 
v = 0.0, a = -1.0 
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R- Function R(0.75,v t - 1,0.1) R-F unction R<0. 25, v,- 1,0,1) 
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R-Function R(0.75,0,a,0,t) 



Figure 3c. Effect of a on R 015 0 (a,0,t) 
q = 0.75, v = 0.0 


Figure 3d. Effect of a on R LOO O (a,0,t) 
# = 1.00, v = 0.0 


NASA/TP— 1999-209424 


6 



Eigen-property 

The R -function also has the eigenfunction character under qth order differintegration with 
v = 0. This is seen as follows. Consider 






„=o r((n + i)?) 

Now, Oldham and Spanier (1974 p.67) prove the following useful form 


’ A 1 r(p-v+i) 


p>-\. 


Applying this equation we have 

d q R (act) = yMfcfJ 

^ rfo) 


nq-\ 


q > 0. 


71=0 


Now let ft = m + 1, then, 


=i 


or 


r((m + l)?) 

(a) m (r-c) (ra+0? - 1 


9 > 0 


C d^ R q ^{ci,cj s ) {a)R q 0 (cz,c^)h- a Hni 


q>0. 


[m + \)q) 

The right most term in equation (26) is zero for t ^ c , thus, for t > c the final result is 

c d q R v0 {a,c,t)= aR q0 {a,c,t) t > c, q> 0. 

Thus, for a = 1 the function is seen to return itself under qth order differentiation. 
Differintegration of the /?-Function 

It is of interest to determine the differintegral of the R -function, that is 

. K ,X r ((„ + i),-v) § r(( n + i>,-v) 

Oldham and Spanier (1974 p.67) prove the following useful form (equation (23) repeated) 

r -jp r(p + l)[x-a] /, " v 

8 x * -fl “ r(p-v+ 1) p> ’ 

which is applied to equation (28) to yield 

°° (n\ n (t 

7 — 4 - v > o. 


So r ((« + 1 )?-( v + w )) 


Thus we have the useful result 


c d“R, v [a,c,t]= R qXv+u) [a,c,t] 


q>v. 


That is, u order differintegration of the R -function returns another R -function. 
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Relationship Between R q m and R q 0 
From the definition of R we can write 


^0 r{(n + \)q-mq) 


r{{n-m + \)q)' 


Letting n-m = r, yields 




r((r+i)?) 


(32) 


(33) 


or 


*"b-'->)-<pr± ^ r X:'?T + w i (a),( ' ' c)l " ,h " 


£ r((r + l)9) ' r((r + l)?) 

Recognizing the first summation on the right hand side as R q0 (a,c,t ), gives the final result as; 


(34) 


R q . mq (a,cj)= (a)"' R fJt (a, c,t)+ (a)"' £ ~ • 

r=-m I \v 1/9/ 

It is noted, that when (r+ l)<? < 0 and integer the elements of the summation term vanish. 

Fractional Impulse Function 

Consider the function R 0 (0,0,r), then we can write, 

«*» / y* +( rt+ 

R q .o (0.0, t) = lim R 0 (a,0, t) = lim X 777 ppr- 

a ^° M 0 Sr((n + l)?) 

In the limit the terms n > 0, of the summation vanish, thus 

( n 'ft q ~ { t q ~ x 

From equation (19) the associated Laplace transform pair is given by; 


(35) 


^K.o(0,0,r)}=- 


Re(g)>0, Re(s)>0. 


(36) 


(37) 


(38) 


Relationship of the /?-f unction to the Elementary Functions 

Many of the elementary functions are special cases of the R -function. Some of these are 
illustrated here. 


Exponential Function 

Consider R lo (a,0,t), by definition, we have 


r(n + 1) ,, = o n! 


thus 


R uo (aAt)=e a '. 


(39) 

(40) 
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Sine Function 

Consider aR 2 Q {- a 2 by definition, we have 


aR 


~ ( 2 V .(rt+l)2-I r 2 .2 4.5 

/ 5 ~ \ X- l-a J r ' a 1 a 1 

, n (-a',0,f)=a> - — t 1 7 = a{t 1 .. 

• oV ' £ (2« + l)! 1 3! 5! 


thus 


a R 10 {-a 2 ,0,r)= sin (a/). 


Cosine Function 

The cosine function relates to R 2] (— a 1 again by definition 

- ( V*(n+I)2-1-1 ~ / „2Y,2* f 2.2 4,4 

R 7l (-a 2 ,0,t )=ft£lL r = i ( ~r V = 1- — + — 

2jV ’ ' to r((n + 1)2-1) h {In). \ 2! 4! 


thus 


i? 2 i(-a 2 ,0,r)= cos(ar). 


Hyperbolic Sine and Cosine 

Consider aR 2 0 (# 2 ,0,r), by definition, we have 

**“*■' A ' )- - 1 


_ 3.3 5.5 

at at 


~ \at 1 h.. >, 

3! 5! I 


thus, 

In similar manner 

/^-Function Identities 


aR 2 o (a 2 ,0, t)=s\nh(at). 


R 2l (a 2 ,0,t)= cosh(ar). 


(41) 

(42) 


(43) 

(44) 


(45) 

(46) 

(47) 


Trigonometric Based Identities 

A number of identities involving the R -function may be readily shown based on the 
elementary functions. The exponential function, equation (40) 



R l0 (a, 0,x) = e ax , 

(48) 

may be expressed as 

H 

II 

>3 

o 

© 

* 

(49) 

Then from equation (42) 

sin (ax) = a R 10 (- a 2 ,0, x) 

(50) 

and expressing the sine function in complex exponential terms gives 



sm(x) = -(e‘ x -e- ix ). 
w 2 1 ’ 

( 51 ) 
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Combining equations (49) ,(50) and (51) then yields the identity 



R, 0 (- 1,0, x) = L (R, 0 (/,0,x) -R, o (- /,0, x)) 
2/ 

(52) 

In similar manner using the cosine function, equation (44) 

cos(x) = R 2i (- 1,0, x) = I (e ix + e ' ix ) 

from which 

(53) 


| 1 

R 21 (- 1,0, x) = -(/?, 0 (z,0,x)+ R l0 (- z,0, x)). 

M 

(54) 

The hyperbolic functions may also be used as a basis, using sinh function, yields 


/? 20 (l,0,x)= l(/?, 0 (l,0,x)- R ]0 {- 1,0, x)) 

(55) 

The cosh function gives 



^i(l,0,x)=I(/? Ii0 (l,0,x)+/? 10 (-l,0,x)) 

(56) 


Many other identities may be found based on the known trigonometric identities, a few examples 
follow, from 

sin 2 (x)+cos 2 (x)= 1, 


we have 


From the identity 


derives 


Rio (- 1,0, x)f R ; ; , (- 1,0, x)= 1. 


sin(2x) = 2sin(x)cos(x), 


Ri,o (“ l,0,2x) =2 R 2 0 (- 1,0, x)R 2 , (- l,0,x) 


From the trigonometric identity 

sin(3x) = 3 sin(x)- 4 sin J (x), 

we determine the identity 


R 2 ,o (~ 1,0,3.x) =3R 20 (- 1,0, x)-4/? 20 (- 1,0, x). 


(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


Further Identities 

Other identities may be derived as follows. Let v = q- p, then the Laplace transform of the 
R -function may be written as 

j 


L { R ,. q - P { sq+a s P-i{ s i +a y 


(63) 
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This may be rearranged to give 


Inverse transforming gives the identity 


^(s'+a)' 


! + ■ 


— a 
s q + a 


(64) 


R q . q - P (- a,0,t) =R pfi (0,0, 0- aRj- 




(65) 


Another set of identities follows by factoring the denominator of Laplace transform, thus 




= 5 


1 


-a 

Now a partial fraction expansion of the denominator gives 


(s ql2 -a u2 ls q,2 +a u2 y 


( 66 ) 


= s 


1 

1 


1 , 

^ s' 

2<a 1/2 

2a vz 


2a n 5 

2a ]/1 

s q, 2 -a m 

s q,2 +a 1/2 


s ql2 -a 1/2 

s q/ 2 +a 112 

L J 



_ 


(67) 


Taking the inverse transform, yields 


( 68 ) 


Very many more such identities are possible, indeed because of the generality of the R -function, 
powerful meta-identities may be possible. 



Relationship of the /?-Function to Other Functions 

The generality of the R -function allows it to be related to many other functions. In this 
section it will be related to the important functions discussed in the introductory section of the 
paper. The Laplace transform facilitates determination of the desired relationships. The double 
arrow will be used to indicate the transform pairs, thus for the R -function; 

R (a,c,/)<=>-7 - — Refo-v)>0 (69) 

s —a 


Mittag-Leffler’s Function 

The Mittag-Leffler function and its transform relate to the L-function as; 


E q [-at“]< 


s H + a 


^ R q , q - }( , f) 


The time domain relationship is 

(- «.o. <K [-"’]-£■ 


Also, because [— 


a,c,f] = ^.( v+a )[-a>c,tJ it follows that 
, dr' R q .o ( ~a,c,t)= E q [- a(t - c) q J 


(70) 

(71) 


(72) 
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Argarwal's Function 

The Argarwal function and its transform relate to the R -function as follows; 

? O’ 0 ’/). 




s g - 1 


The time domain relationship is 




nq-\+p 


i^T{nq + p) 


(73) 

(74) 


Erdelvfs Function 

The relationship between the Erdelyi function and the R -function is given by 


(75) 


Robotnov and Hartley Function 

The F-function and its transform relate to the R -function as follows; 

1 

T~~ « 
s +a 

The time series common to these functions is given as; 


F a [~a,t] 


<=> 




_ 00 f 

^(-«A0=^[-^d=L r((n+l)q) ’ 


(76) 


(77) 


n - 0 


Miller and Ross’s Function 

The Miller and Ross function and its transform relate to the R -function as follows 


E, (v,a) <» — — «> R u _ v (, a,0,t ). 


s — a 


The time series common to these functions is given as; 


(a, 0,r)=£,( v ’' a )=S 


0)V 


^r(n + v + l) 


(78) 


(79) 




Example - The Dynamic Thermocouple 

This problem was introduced originally in Lorenzo and Hartley 1998, and frequency domain 
solutions are presented there. Here, it is desired to determine the time domain dynamic response 
of the thermocouple, figure 4, which is designed 
to achieve rapid response. The thermocouple 
consists of two dissimilar metals with a common 
junction point To achieve a high level of 
dynamic response, the mass of the junction and 
the diameter of the wire are minimized. Because 
the wires are long and insulated they will be 
treated as semi-infinite (heat) conductors. This 
analysis will determine the time response of the 
junction temperature T h (s) in response to the 


Q(f) 


i 


■k 2 ,a 2 


Q 2 


t) 


Figure 4. Dynamic Thermocouple 
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free stream gas temperature T {s) . For the semi-infinite conductors the conducted heat rate 
Q{t) is given by 

Qk)=-r=A n T b , ( 8 °) 

V°7 

where k is the thermal conductivity and a is the thermal diffusivity. For the transfer function 
the effects of initialization are not required, therefore, all y/{t)' s are zero. Thus the following 
equations describe the time domain behavior: 

Q i (t) = hA{T(t)-T b (t)} (81) 


T b (t ) = — o d;' (a {t)-QAt)~Qi (0) 

WC 


(82) 


G, (' ) = -== o D) ' 1 T, (< h ~r= (o 4'' !7 s (' ) + V, . X , o .0. 1 )) and (83) 


e ! (,)=-^= t D, l ' ! r,( I )=4i=U, l ' ! r 1 (,)+r J (j;.x.->. «.»)) (84) 

Ja? Jet, 


where h A is the product of the convection heat transfer coefficient and the surface area and 
wc v is the product of the junction mass and the specific heat of the material. Taking the Laplace 
transform of these equations yields 


_1 

wc 


V \ ' 


Q,(s) = hA{T e (s)-T b {sj) 

(85) 

~ \Q; (s)~ Q { (*)- Q 2 (5)]+ Wi (0 

\ 

/ 

(86) 

<2,(0 = A -(? u ‘T b {s)+Vi(s)] 
V«i 

1 

(87) 

Q 2 ( s )= rM* T b {s)+yf 2 ( s l 

V«2 

> 

(88) 


Eliminating the Q's, and solving for T b (s) yields 


T b (s)= 


f A 

WCy 

s+bs y 2 + c 


hAT, (s) — ( s ) — //,(s)+ sy,(s) 


: (89) 


where b = 


wc,. 


Ic 1 c 

*1 { k 2 


a, jet. 


hA 

, and c = . Factoring the leading denominator and 

wc 
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expanding in partial fractions gives 




f _L 


I '\ 


U’C. 


P: ~ A | Pi~Pz 


5 ,/2 +A 5 1/2 + /3, 


hAT(s) — S=y/-,(0 — ^=l/A 2 (5)+5V/ 3 (s)k (90) 


where A = — + -^4b>^-Ac and A = — — — v/A -4c . Then with appropriate choices for the 

functions of s in the right most bracket this equation may be inverse transformed to yield the time 
domain response. To demonstrate the value of the R -function, we select (determine) 

27’fO) 1 

y/ 3 (s) = T b (0 )/ s, Further assume T g (f) = 2 T b (0) + 1 => T g (s) = — - — - + — , and 
\j/ l (.v) = y/ 2 (.v) are arbitrary functions of time. The solution may be written directly as: 


b 1 


T b (t) = 


hA 


[ K (0)/? 1/2 .. ( (- p x ,o ,t)+ r U2 '_ 2 (- A ,0,0 


wc„03 2 -A) 

- 2T b (0 )/?,*., (- p 2 ,0 1,0- A/ 2,2 (- Pi A 0 ] 

J"{^i/2,o( — Pi A? — f) — A/ 2,0 ( — Pi >0, t — t))// , {r^dx 


wcAPi -Pi) 


k k 
*1 . . 2 




+ 


W C v (p2~Pl) 


T b (ofe/2.0 (- Pi A 0 - A/ 2 , 0 (- Pi AO} 


(91) 


Further Generalized Functions 

Functions yet more general than the /? -function may be developed. One such function will 
be derived here. It is simpler here to work backward from the s-domain to the time domain. Thus, 
we consider the following function 


G(0 = 




where and rare not constrained to be integers. Then this may be written as 

Y r 


G(s) = s'” 


qr 


f i-± 


) 


Now the parenthetical expression may be expanded using the binomial theorem to give 


g«=* 


r(l-r) f-a 
fer(l + j)T(l-j-r)l s' 


V 


<l. 


or 


G(s)= I 


r(l-r) 


^r(i + y)r(i -;-0 


(-ays 


v-qr-qj 


(92) 

(93) 

(94) 

(95) 
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This expression may be term by term inverse transformed yielding 


g, „ ,m= z 7 r w r)( t ■ Re( «'' ■ ^ a Re(s ^ °- 4 < 1 ■ 

Thus we have the following transform pair 

4{G„,v,rk f ]}= / ^ y ' Re{qr-v)>0, Re(s)>0, >0. (97) 

(s*-aj 5 

The form of equation (96) presents evaluation difficulties, since when r is an integer 
r(l — r) and F(l — j — r) can become infinite. Equation (96) maybe rewritten as follows: from 
Spanier and Oldham (p.414, eq.43:5:5) 

r (x-n)=? ~^ \~7 r B = 0 ’ 1 ’ 2 ’- (98) 

(x-lXx-2)---(x-n) 

where (l - x) n is the Pochhammer polynomial. From this result with x = 1 - r , we can write 

Itj-j-r)-, v -MSbd y = l,2, (99) 

(-rX-l-r)'O-r-j) (r), 


Substituting this result in equation (96), yields the following more computable results 
~ r 1 y {(- r\- 1 - r>- ■(!■ - T -OX-y 

*'' 1 ’ 4 r ( i + itdr+ ib-v) 




or in terms of the Pochhammer polynomial 

[ - ' ] = ■ Re ( - - 0 " °. R e ( , } ^ 0 . | ^ | ‘ . (,00 


In similar manner relationships of increasing generality may be determined. Podlubny (1999) 
presents a form that is a special case of the G -function where r is constrained to be an integer. It 
is also clear that taking r = 1 specializes the G -function into the R -function. It is the authors’ 
judgment that the F- and R -functions will prove to be the most useful in practical applications. 
Table 1 summarizes the advanced functions studied in this paper along with their defining series 
and Laplace Transforms. 
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Table 1 Summary of Defining Series and Laplace 


Function 

Time Expression 

Laplace Transfom 

Remarks 

Mittag-Leffler 


/ 

5 ( 5 ^ - a) 

0 differintegral of 

Agarwal 

Ea ^ ^ l?r (am+p) 

s a -l 


Erdelyi 

°° jtn 

^a,0 0) — j-,/ n\ 

^r(am+j8) 

v> Tfm+i) 

a,p >0 

Robotnov / Hartley 

**£*♦<» 

1 

s q -a 

eigenfunction 

Miller-Ross 

- * *+V 

£,(v,a)=]T , 

s~ v 

s-a 


Current Paper 

- /i, (tt-H L-l-i’ 

s v 

s q -a 

eigenfunction & differintegral 

Current Paper 

i-o (- 0 ' r(l + /X((r + j)q - v) 

V 

5 

(s’ -aj 

eigenfunction Sc differintegrals 
repeated & partially rep. 


Summary 

This paper has presented a new function for use in the fractional calculus, it is called the 
R -function. The R -function is unique in that it contains all of the derivatives and integrals of 
the F-function. The R -function has the eigen-property, that is it returns itself on q\h order differ- 
integration. Special cases of the R -function also include the exponential function, the sine, 
cosine, hyperbolic sine and hyperbolic cosine functions. Further, the R -function contains, as 
special cases; the Mittag-Leffler function, AgarwaFs function, Erdelyi's function, Hartley's 
F-function, and Miller and Ross's function. Numerous identities are possible with the 
R -function some of these have been shown in the text. 

The value of the R -function is clearly demonstrated in the dynamic thermocouple problem 
where it enables the analyst to directly inverse transform the Laplace domain solution, 
(operational (s) form) to obtain the time domain solution. 

A further generalization of the R -function, called the G -function brings in the effects of 
repeated and partially repeated fractional poles. This generalization carries increased time 
domain complexity. 

A R -function based trigonometry is also possible. It is a generalization of the conventional 
trigonometry, and will be the subject of a future paper. 
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